Modulational and Parametric Instabilities of the Discrete Nonlinear Schrodinger 
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We examine the modulational and parametric instabilities arising in a non-autonomous, dis- 
£NJ , crete nonlinear Schrodinger equation setting. The principal motivation for our study stems from 

5^ ' the dynamics of Bose-Einstein condensates trapped in a deep optical lattice. We find that under 

C^H, periodic variations of the heights of the interwell barriers (or equivalently of the scattering length), 

additionally to the modulational instability, a window of parametric instability becomes available to 
the system. We explore this instability through multiple-scale analysis and identify it numerically. 
Its principal dynamical characteristic is that, typically, it develops over much larger times than the 
modulational instability, a feature that is qualitatively justified by comparison of the corresponding 
instability growth rates. 



I. INTRODUCTION 



The modulational instability (MI) is a general feature of discrete as well as continuum nonlinear wave equations. For 
this instability, a specific range of wavenumbers of plane wave profiles of the form u(x, t) ~ exp(i(kx — tot)) becomes 
unstable to modulations, leading to an exponential growth of the unstable modes and eventually to derealization 
q . (upon excitation of such wavenumbers) in momentum space. That is equivalent to localization in position space, and 
' hence the formation of localized, coherent solitary wave structures [1]. 

The realizations of this instability span a diverse set of disciplines ranging from fluid dynamics [2] (where it is 
usually referred to as the Benjamin- Feir instability) and nonlinear optics [3] to plasma physics [4]. One of the earliest 
contexts in which its significance was appreciated was the linear stability analysis of deep water waves. It was much 
later recognized that the conditions for MI would be significantly modified for discrete settings relevant to, for instance, 
the local denaturation of DNA [5] or coupled arrays of optical waveguides [6,7]. In the latter case, a relevant model 
. is the discrete nonlinear Schrodinger equation (DNLS), and its MI conditions were discussed in [8]. Most recently, 
the MI has been recognized as responsible for dephasing and localization phenomena in the context of Bose-Einstein 
condensates (BEC) in the presence of an optical lattice i.e., a sinusoidal external potential [9-12]. 

In the context of BECs which are among the principal motivations of this work, another interesting possibility arises. 
For a "deep" optical lattice (i.e., if the wells of the spatially periodic potential are well-separated and sufficiently 
high), it has been shown that the relevant mean field model that describes the behavior of the condensate, at 
T = 0, is the discrete nonlinear Schrodinger (DNLS) equation [10,13-15]. The optical lattice can be created by two 
I ■ counterpropagating laser beames forming a standing wave interference pattern. 

Our interest in the present work is in introducing an explicit temporally periodic modulation in the coefficients of 
the DNLS and examining the instabilities that may arise (for uniform solutions) . In the BEC setting, there is a number 
of potential realizations of such a non-autonomous DNLS equation. For instance, the heights of the interwell barriers 
of the optical lattice are proportional to the intensity of the lasers, and can be easily periodically modulated in time. 
This induces an oscillating tunneling amplitude of the condensates between adjacent wells, as well as an oscillating 
interaction energy of the condensates trapped in each well. An alternative possibility involves the periodic modulation 
of the scattering length of the interaction between the atoms via a Feshbach resonance, i.e., an external magnetic 
field; see e.g., [16]. The possibility of this, so-called, Feshbach resonance management (FRM) of the interaction has 
generated a large interest recently due its robust effect on patterns, coherent structures and its potential for avoiding 
collapse, see e.g., [22,17-20]. 

In this short communication, we revisit the modulational instability criteria in the DNLS equation (which were 
originally derived in [8]), but in the presence of the periodic modulation of the DNLS tunneling and interaction 
parameters, as motivated above. We choose the simplest possible periodic modulation (a sinusoidal variation of 
the atomic scattering length) and derive the modulational stability equation, which in this case becomes a modified 
Mathieu equation. In the absence of the periodic perturbation we recover the results of [8]. In the presence of 
such a term, an additional, parametric instability becomes possible. These new domains of instability appear due 
to parametric resonance, whenever the parameters of a system vary periodically with time. In contrast to ordinary 
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resonance, where we have growth proportional to the time variable, here the growth is exponential, as in the case of 
customary modulational instability; see e.g., [21]. We implement a multiple-scale expansion to identify the instability 
domain boundaries and subsequently numerically examine our analytical predictions. The numerical investigations 
indicate that this parametric resonance sets in over much longer time scales than the modulational instability (for 
the same perturbation amplitude). Similar considerations for the continuum problem can be found in [22,23] .In [22] 
a slightly more restricted ansatz was used [q = was used in our Eq. (5) below], and in [23] it was the dispersion 
that was varying with time, instead of the nonlinearity strength as in our case. The continuum limit of our analytical 
results (modulo the relevant rescalings/adjustments) has been found to agree with the results of [22,23]. 

Our presentation is structured as follows. In section II we present the mathematical framework and our analytical 
results. In section III we corroborate these results with numerical simulations. Finally, in section IV, we summarize 
our results and present our conclusions. 



II. SETUP AND ANALYTICAL CONSIDERATIONS 



We study the discrete nonlinear Schrodingcr equation in the form: 

ij> n = -D(t){tp n+1 + - 2V>„) + a(i)|V„| 2 Vn + E(t)ip n , (1) 

where the coefficients D(t), a(t) and E{t) arc time dependent. If one sets ip n (t) = <f) n {t)e~ l §o ( 2D ( S )+ E ( s )) ds ^ ^hen ^ 
is reduced to 

i]> n = -D{t){4> n+1 + <f) n -i) + a(t)\(j) n \ 2 (b n . (2) 

Now, if we additionally set r = r(t), </>„(t) = 4> n (t), D(t) = D(t), a(r) = a(t) and choose ^ = D(t), then (2) is 
equivalent to 



i(f>n = -(4>n+l + <j>n-l) + b(t)\(j) n \ 2 4> n , (3) 



where b(t) — jjth, and we must notice at this point that, with a small abuse of notation, we will continue to use t 
instead of r. From now on, we will study (3), where we assume that 

b(t) = l + esm((jt), (4) 

as motivated earlier. It is easily verified that 

/.x i(2cosqt—f b(s)ds+qn) /r x 

v n {t) = e Jo (5) 

where q is the wavenumber, is an exact solution of (3). In order to examine the modulational stability of this plane 
wave solution, we use the ansatz 

Mt) = Vn(t) [1 + e {a(ty kn + (3{t)^ lkn )] , (6) 

where k is the perturbation wavenumber and a(t), (3(t) are complex, time dependent fields. Substituting (6) into (3) 
and keeping only the 0(e) terms, we obtain the following first order, coupled system for a and /?*, where * denotes 
the complex conjugate: 

J a = 4sin(^-^)sin(^)a + 6(a + /3*) (7) 
0* = 4 S in(^-^)sin(^)/r -b(a + p). (8) 
These can be combined to give a second order equation: 

a + (2iB - h -)a + (A 2 - B 2 + 2Ab - i^(A + B))a = 0, (9) 

where 
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A = 2cosg(l - cosfc) (10) 
B = 2 sing sin fc (11) 

are constants that depend only on the wavenumber q and the perurbation wavenumber k. It is worth noticing that if 
e = in (4), then the modulational instability criterion reads A(A + 2) < which is in accordance with Eq. (10) of 
[8]. Next, we substitute (4) into (9) to obtain, up to order 0(e), the equation: 

&+(2iB - euj cos(ut))a + (A 2 - B 2 + 2A + e(2Asm(ut) - iuj(A + B) cos(wi))a = 0. (12) 

This is an equation with periodic coefficients [24] , hence it is natural to examine whether parametric instabilities may 
arise due to the temporal modulation. In view of this perturbation, if we use a regular perurbation expansion of a as 

a series in e, a = a a + eat\ + 0(e 2 ), we find that secular terms exist if A = Aq = — 1 + \Jl + which is equivalent to 

q = qo = arccos( 2 ^ A " osk ^ )■ Thus, it is necessary to implement a multiple scale analysis [25], expanding q as a series 
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q = qo + e qi +0(e z ). (13) 
After lengthy but straightforward calculations, it is found that the value of q\ is 

?i = ± - V " • (14) 

4^1 + ^(1 -cosfc) sin q 

After the substitution of (14) into (13) and using the expression (10) for A, we obtain the boundaries of the instability 
domain on the (q, e) plane to be: 



4W1+ ^-(1 - cosfc) sin q 
±e=-Y 4 - == (q-qo). (15) 



III. NUMERICAL RESULTS 



In order to check the validity of our analytical approach, we have performed numerical simulations of the equations 
of motion (3) using a fourth order Runge-Kutta scheme. The parameters of the system have been chosen to be 
e = 0.01, u = 1 and fc = ir/2. The initial condition, in accordance with Eq. (6), is a modulated wave 

u n = e lkn + te l{q+k)n . (16) 

The simulations have been performed with a chain of N = 1380 sites, with periodic boundary conditions so that the 
wave numbers q(k) defined modulo 2tt in the lattice, are of the form q = 2-kt/N (fc = 2nR/N), where r (R) is an 
integer (see also [8]). 

Fig. 1 shows the windows of parametric instability represented by the stability limits of Eq. (15). For wavenumbcrs 
between the two straight lines, parametric instability is theoretically predicted from the results of Section II. To 
illustrate this point, we select 4 wavenumbers, one of which (q — 1.50705) is both modulationally and parametrically 
stable; the second (q = 1.51161) is parametrically unstable; the third (q = 1.51616) is again in the window of stability, 
while the fourth one (q = 1.58446) is modulationally unstable. 
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FIG. 1. The plane of (e, q) where the region of the parametric instability represented by Eq. (15) is shown by two straight 
lines. Wavenumbers between the two "stability limits" will be (to leading order) parametrically unstable. We also highlight 4 
wavenumbers for which we conduct numerical simulations, namely q = 1.50705, q = 1.51161, q = 1.51616 and q = 1.58446 , 
that lie on the line e = 0.01 (see text). We use PS (PU) to indicate the parametrically stable (unstable) regions, respectively, 
and MU to indicate modulationally unstable ones. The vertical line q — tt/2 is the boundary between the modulationally 
stable-unstable regions. 

In the left panel of Fig. 2, we indeed observe at time t = 800 that the perturbed wave is both parametrically and 
modulationally stable. On the other hand, the right panel displays q = 1.51616, which should also be parametrically 
stable according to our linear theory. However, for longer times than the ones displayed in Fig. 2, the latter 
wavenumber has been observed to become unstable. We conjecture that this instability is due to a higher order 
parametric resonance, not captured within the theory of section II. 




FIG. 2. Both panels show the case where we have modulational and parametric stability. The left corresponds to q = 1.50705. 
and t = 800 and the right to q = 1.51616 and t = 800. The squared amplitude (\u\ 2 ) of the solution is shown in the panels as 
a function of the lattice site n. 

The case of Fig. 3 displays the time evolution of the wavenumber predicted to be in the window of parametric 
instability (q = 1.511616). We indeed find that at times t = 570, the configuration becomes parametrically unstable. 




Finally, in Fig. 4 we show the time evolution of a modulationally unstable case for q = 1.58446. One clear feature 
that distinguishes the parametric instability shown in Fig. 3 from the modulational one of Fig. 4 is the time scale 
of the instability development, or equivalently the corresponding instability growth rate. The latter appears to be 
much larger in the case of the modulational instability; as a result, the MI sets in much sooner for equal magnitude 
perturbations of the two plane waves. 
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FIG. 4. The case which is modulationally unstable. The left panel shows the square modulus of the field, before the onset 
of the instability at t — 17, and right after the onset, at t = 56. In this case q = 1.58446. 

From the growth rate predictions of Eq. 10 of [8] and from our analysis, it can be calculated (in the cases shown) 
that the theoretical growth rate in the case of the modulational instability is approximately 0.23218 (forg = 1.58446) 
and the corresponding growth rate in the case of the parametric instability is approximately 0.00115 (for q = 1.51161). 
The considerably larger growth rate of the MI qualitatively justifies the earlier temporal development of the latter. 



IV. CONCLUSIONS 



In this short communication, we have examined the potential of time-dependent coefficients, in the context of a 
non-autonomous discrete nonlinear Schrodinger equation, to induce a parametric instability of plane wave states. 
Such an instability has been identified and its boundaries established at the level of a leading order theory within 
a multiscale expansion. The theoretical findings have been numerically tested through direct simulations and have 
been found to be in agreement with the theoretical predictions (except for the regime between the parametric and 
modulationally unstable wavenumbers, where higher order parametric resonances may ensue). 

The distinctive feature of the present parametric instability is that it has a much longer threshold time for its 
dynamical development in comparison with the modulational instability. This can be both quantified by means of 
the comparison of their respective growth rates, as well as observed in the course of direct numerical simulations. 
Studying the effects of parametric instabilities in other contexts including that of coherent, nonlinear wave structures 
in BEC is very relevant and will be addressed in future publications. 

PGK gratefully acknowledges support from NSF-DMS-0204585 and from the Eppley Foundation for Research. 
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